The magnetic response of disordered metallic rings: the large contribution of the far 

levels 
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We calculate the orbital linear magnetic response of disordered metallic rings to an Aharonov- 
Bohm flux using the BCS model for attractive electron-electron interaction. The contribution of 
all levels including those up to a high energy cutoff results in a much larger value than previously 
obtained using the local interaction model. The possible relevance of our results to the resolution of 
the discrepancy between the experimental and theoretical values for the ensemble-averaged persistent 
currents in these systems is discussed. 
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One of the remarkable phenomena of mesoscopic 
physics is the existence of equilibrium persistent currents 
in small normal metal rings, in the regime where the 
clastic mean free path I is much smaller than the ring's 
circumference L. This was both predicted theoretically 
[1] and observed experimentally [2-7]. While the phe- 
nomenon and the periodicity in the flux are well under- 
stood by now, the magnitude of the ensemble-averaged 
persistent current found experimentally is much larger 
than that obtained theoretically using the model of non- 
interacting electrons [8,9]. An attempt [10,11] to account 
for this discrepancy by the inclusion of electron-electron 
("e-e") interactions [12] indeed increased the theoretical 
value, but still came short by a factor of about 5 [13]. 

Along with the large value of the ensemble-averaged 
persistent current observed, its sign poses significant 
questions as well. While in the first experiment [2] 
the sign of the ensemble-averaged persistent current was 
tentatively identified as diamagnetic, it was not deter- 
mined definitively (see Ref. [13]). Later experiments re- 
ported a predominantly diamagnetic ensemble-averaged 
persistent current [5-7]. Theoretically, the paramag- 
netic/diamagnetic sign follows for repulsive/attractive ef- 
fective e-e interactions (the interplay between the renor- 
malized [14] repulsive e-e interaction and the phonon me- 
diated attractive interaction determines the sign of the 
effective interaction [13]). Thus, the diamagnetic sign of 
the persistent currents observed in Gold [6], Copper [2], 
and Silver [7] suggests that these materials are very weak 
superconductors [15]. 

The model used to treat both signs of the interaction 
was that of local e-e interaction. The discrepancy be- 
tween the theoretical values and the experimental results 
exists for repulsive [10], as well as attractive interactions 
[13]. For the repulsive case, the higher order renormal- 
ization which reduces the first order result is responsible 
for this discrepancy. For the attractive interaction case, 
the higher order renormalization does not reduce, and 
even somewhat increases the first order result. However, 
the restriction on the bare interaction constant, due to 
the experimental fact that Gold, Copper and Silver are 
not superconducting in observed temperatures, limits the 



value of the theoretical result in this case. For the mag- 
netic response (the derivative of the persistent current 
at zero flux), there exist a further logarithm in the nu- 
merator [10], irrespective of the above renormalization 
of the effective interaction. Here we show that the use 
of the BCS Hamiltonian to model attractive e-e interac- 
tion gives rise to a much larger cutoff of this logarithm, 
and hence to a significantly larger value for the mag- 
netic response to flux. The enhancement of the magnetic 
response is a result of the contribution of high energy 
levels (to be called "far levels" ) originating from pairing 
correlations which persist up to the energy of the De- 
bye frequency wd. Thus, correlations on an energy scale 
much larger than the Thouless energy affect the magnetic 
response. 

The contribution of the far levels was investigated 
in connection to superconductivity in small grains [16]. 
There, the correlations of levels much further than the 
superconducting gap from the Fermi energy E-p, up to 
cl>d prove to be significant. Using an exact solution of 
the reduced BCS Hamiltonian [17,18], it was shown that 
this contribution results in a much larger condensation 
energy than that given by the BCS value, in a wide pa- 
rameter regime in which superconducting correlations are 
well developed. The far levels also affect single particle 
properties of a superconducting grain [19], such as the 
Matveev-Larkin parameter [20] . 

The difference between the local interaction model and 
the BCS Hamiltonian lies in the q dependence of the in- 
teraction in each, where q is the sum of the incoming 
(outgoing) momenta of the pair scattering. While the lo- 
cal interaction model assumes that the bare interaction 
does not depend on any momenta, and specifically on 
q, the BCS Hamiltonian assumes a sharp cutoff at small 
total momentum, and considers the interaction to be a 
delta function at q = 0. While for physical properties 
related to the superconducting phase of the system, near 
and below T c (such as T c itself) the two models give sim- 
ilar results, we show that for the magnetic response in 
the perturbative regime (equivalent to T > T c ) the re- 
sults in the two models differ substantially. We discuss 
the relation between the two models, and the effect of 
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relaxing the BCS assumption by taking the interaction 
to have finite width in q. 

Recently it was suggested that the largeness of the ob- 
served current is due to AC noise [21,22] or interactions 
of the electrons with two level systems [23], and a rela- 
tion to dcphasing was suggested (one should also keep 
in mind that the latter interactions may lead to an ad- 
ditional attractive interaction). We, however, consider 
T = 0, and do not include dephasing (we consider the de- 
phasing length Lff, to be larger than all relevant lengths) . 
A finite will result in a suppression of the persistent 
current, exponential in L/L<p. 

We consider a quasi one dimensional disordered ring 
penetrated by a constant Aharonov-Bohm (AB) flux $ in 
its middle. We calculate the magnetic response, denoted 
by i<2, to first order in the e-e interaction, for attractive 
interaction using the reduced BCS Hamiltonian. 

The result we obtain with logarithmic accuracy is 
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where A is the dimcnsionlcss pairing parameter [see 
Eq. (6)], d is the level spacing, i?Th = 1 2 /(stL 2 ) is the 
Thouless energy, which is the inverse of the time to cir- 
culate the ring (s — 1,2,3 is the effective dimension of 
the ring for diffusive motion), <f> = hc/(2e) and (• • •) 
denotes ensemble averaging. The upper logarithmic cut- 
off E co = min [uip, , 1 /r] represents the fact that energies 
(measured from the Ep) further than £/rh, up to E co (far 
levels) contribute to the magnetic response (r is the elas- 
tic mean free time). Their contribution enhances F2 by 
about an order of magnitude due to the much larger log- 
arithmic cutoff as compared to £>rh of the known result 
[10] [for 1/(t-Et1i) ~ 10 4 the ratio of the logarithms is 
roughly 10]. 

In order to obtain Eq. (1) we consider the general 
Hamiltonian 
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where U(r) is the external potential which includes the 
disorder, and H cc represents the e-e interaction. In 
the London gauge A = <J>/(27r /?)</> where p is the dis- 
tance from the origin and <fi is in the clockwise direc- 
tion of the ring. The free energy of the system and 
the persistent current are flux dependent, and related by 
/($) = -dF/d$. By time reversal symmetry 1(0) = 0, 
and for small flux 
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We now turn to the calculation of F 2 to first order 
in the interaction. Pcrturbative analysis of the reduced 
BCS attractive interaction [Eq. (6) below] is valid, at 



T = 0, for A < l/\n(uip,/d) [16] and at finite tempera- 
ture (T > d) for T > T c . We take as the unperturbed 
Hamiltonian 
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and the magnetic field together with the e-e interaction 
as perturbation 



Hi — H ee + 
E / dr ^{r) 
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We calculate perturbatively the energy to third order in 
Hi, and consider only the terms which are second order 
in the flux and first order in the interaction. 

We denote by \i) the eigenstates of the noninteracting 
electrons in the disordered ring without magnetic field. 
In this basis, within the reduced BCS model, we obtain 
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Here Cj destroys an electron in the state |i) with wave- 
function Xi{ r ) an( l Pij = is the matrix element 
of the momentum parallel to the ring's direction. The 
second sum is restricted to levels within lo-q of Ep. Note 
that contrary to the usual convention, A is negative for 
attractive interaction. We assume that the width of the 
ring is much smaller than its radius [24] . The A term is 
then interaction independent. We choose the Xi' s to be 
real, and then P^ is pure imaginary and Pa = 0. 

Using third order perturbation theory [25] we find that 
for T = to first order in A 
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where < i < wp, represents states whose energy is be- 



tween E-p and Ep + up, , and Wjj = e j — 
disorder averaging over \Pij\ 2 we obtain 
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Performing 



(8) 



which is roughly constant for u>ij < 1/r and zero for 
LUij > 1/r. In Eq. (7), for energies smaller than up, and 
1 /t there is a double sum and a second power of energy 
in the denominator. The result is therefore a logarithm 
with an upper cutoff which is given by min \u>p>, 1/r]. We 
evaluate the sum in Eq. (7) taking for simplicity the non- 
interacting spectrum to be equally spaced, and obtain 
Eq. (1). 

This equation differs from the result of Ambegaokar 
and Eckern (AE) [10] by the upper cutoff of the log- 
arithm, being min [up, , 1 /r] in comparison to £rh in 
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Rcf. [10]. We now explain the origin of this differ- 
ence in some detail. AE consider the same Hamilto- 
nian as in our Eq. (2), only with local e-e interaction 
V(r\ — r-2.) = \N(0)~ 1 5(ri — r^)- In terms of the nonin- 
teracting eigenstates H ee is then given by 

#c°c Cal = \Y. ^ c UW^ ct , (9) 

ijkl.a 

where V ijk i = AiV(O)" 1 J Xt(r)Xj{r)Xk(r)xi(r)dr and 
a ^ a due to the interaction being local. The dependence 
of the energy on flux, to first order in the interaction, can 
be written as 

E V in^)= E +£*««(*)■ 

i,j<0 i^j<V i<0 

(10) 

We now show that the second, diagonal part contribution 
to the magnetic response corresponds to the BCS result 
as given in Eq. (1). AF can be expanded to second order 
in the flux to obtain Fi for the local interaction model by 
following the same procedure done above for the reduced 
BCS Hamiltonian, with H ee given by Eq. (9). The result 
is 

«=*(^) 2 ( E ^Eflto -"«.))■ 

\3<0,i>0 % 3 2<0 / 

(11) 

The diagonal contribution to this result is given by the 
/ = i terms, and since to leading order Van — V{jji = 2Xd 
[26] we obtain for the diagonal contribution in the local 
interaction model the same expression as Eq. (7), without 
the explicit cutoff at wd- Due to Eq. (8), the diagonal 
contribution gives Eq. (1) with E co = 1/r. Therefore, 
the BCS approximation is equivalent to assuming that 
the flux dependence of the offdiagonal matrix elements 
is small, and can be neglected. In the case of the local 
interaction model this is not the case, and the contribu- 
tion of the offdiagonal elements to the magnetic response 
is significant, and opposite in sign to that of the diago- 
nal element. This results in a partial cancellation, and 
therefore a reduction of the high logarithmic cutoff in 
Eq. (1). 

Another way to understand the relation between the 
two models is to consider the q dependence of the in- 
teraction in both (where q is the sum of the incoming 
momenta). AE derive their result as a sum over all q's 
(see Eqs. [(12)-(14)] of Ref. [10]), and the result within 
the BCS Hamiltonian [Eq. (1)] corresponds to their q = 
term. For the local interaction model, the contribution 
to the magnetic response of the high energies, above -&rh 
exactly cancels between equal magnitude and opposite 
signs of the q = term and the sum of all q ^ terms. 
This makes the q independent assumption for the bare in- 
teraction crucial. The existence of excess interaction at 



small total momentum q would thus significantly affect 
the result. The BCS interaction assumes just that, the 
existence of excess interaction at q = 0. 

A physical justification for taking a q dependent inter- 
action can be obtained, for the attractive phonon medi- 
ated interaction, from the usual restriction that all (in- 
coming and outgoing) scattered states are within cjd of 
E-p (see e.g. Ref. [27]). This restriction implies a signifi- 
cant q dependence on the scale of q c = uj-q/vf. One can 
therefore take a q dependent attractive interaction in the 
form of a step function with width q c . For q c < 2tt/L 
one finds that only the q = term survives, and the BCS 
result for the magnetic response is recovered. 

For metals such as Gold, Copper and Silver, this con- 
dition requires the circumference of the ring to be smaller 
than 0.2 — OA/im, which is an order of magnitude smaller 
than the relevant experimental lengths [2,5,7]. The ef- 
fects of relaxing this condition, i.e. having q c > 2ir/L, as 
well as the effect of the q dependence on the persistent 
current itself will be considered elsewhere [28] . 

The central result of this Letter is the large interaction 
correction to the derivative of the persistent current at 
zero flux, within the BCS model, as is given in Eq. (1). 
Physically, this is due to the large contribution of the far 
levels, up to min [wd, 1/t] from the Fermi level. 

Our calculations were done for T = 0. At finite tem- 
perature T < min [wd, (but neglecting dephasing, 
assuming L$ » L) the magnetic response would be given 
by Eq. (1) with d replaced by T. Therefore, on top of the 
large magnetic response, we predict a weak, logarithmic, 
temperature dependence of the magnetic response up to 
T « min [loy> , 1/t] . The weak temperature dependence is 
due to the addition of the many small contributions, all 
having the same sign, of the levels up to the large energy 
cutoff (see also Refs. [16,29]). 

The local interaction model and the reduced BCS 
Hamiltonian are two models used to describe the effec- 
tive e-e interaction resulting from the retarded electron- 
phonon interaction. Though different, these two models 
give similar results for many properties of the supercon- 
ductor near and below T c . The robustness of the physics 
near and below the transition temperature makes the dif- 
ferences between the two models irrelevant. However, the 
magnetic response in the perturbative regime is a more 
subtle property, that distinguishes between the two mod- 
els. 

Our result for the reduced BCS Hamiltonian also fol- 
lows by expanding Richardson's exact solution [17]. This 
has the analytical merit of being a first order expansion 
of an exact solution, as well as an additional viewpoint 
on the contribution of the far levels, and will be given 
elsewhere [30]. 
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